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ABSTRACT 

We present a description of the data reduction and mapmaking pipeline used for the 2008 observing season 
of the Atacama Cosmology Telescope (ACT). The data presented here at 148 GHz represent 12% of the 90 TB 
collected by ACT from 2007 to 2010. In 2008 we observed for 136 days, producing a total of 1423 h of data 
(11 TB for the 148 GHz band only), with a daily average of 10.5 h of observation. From these, 1085 h were 
devoted to a 850 deg 2 stripe (1 1 .2 h by 9? 1) centered on a declination of -52?7, while 175 h were devoted to a 
280 deg 2 stripe (4.5 h by 4P8) centered at the celestial equator. We discuss sources of statistical and systematic 
noise, calibration, telescope pointing, and data selection. Out of 1260 survey hours and 1024 detectors per 
array, 816h and 593 effective detectors remain after data selection for this frequency band, yielding a 38% 
survey efficiency. The total sensitivity in 2008, determined from the noise level between 5 Hz and 20 Hz in the 
time-ordered data stream (TOD), is 32/xK v / s in CMB units. Atmospheric brightness fluctuations constitute the 
main contaminant in the data and dominate the detector noise covariance at low frequencies in the TOD. The 
maps were made by solving the least-squares problem using the Preconditioned Conjugate Gradient method, 
incorporating the details of the detector and noise correlations. Cross-correlation with WMAP sky maps, 
as well as analysis from simulations, reveal that our maps are unbiased at multipoles t > 300. This paper 
accompanies the public release of the 148 GHz southern stripe maps from 2008. The techniques described here 
will be applied to future maps and data releases. 
Subject headings: Microwave Telescopes, CMB Observations 



Departamento de Astronomi'a y Astrofisica, Facultad de Ffsica, Ponti- 
ficia Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile 

2 Department of Physics and Astronomy, University of British 
Columbia, Vancouver, BC V6T 1Z4, Canada 

3 Department of Astrophysical Sciences, Peyton Hall, Princeton Univer- 
sity, Princeton, NJ 08544, USA 

4 Department of Physics and Astronomy, The Johns Hopkins University, 
3400 N. Charles St., Baltimore, MD 21218-2686, USA 

5 Joseph Henry Laboratories of Physics, Jadwin Hall, Princeton Univer- 
sity, Princeton, NJ 08544, USA 

6 Canadian Institute for Theoretical Astrophysics, University of 
Toronto, Toronto, ON M5S 3H8, Canada 

7 Department of Physics and Astronomy, Rutgers, The State University 
of New Jersey, Piscataway, NJ 08854-8019, USA 

8 Department of Astrophysics, Oxford University, Oxford OX1 3RH, 
UK 

9 School of Physics and Astronomy, Cardiff University, The Parade, 
Cardiff, Wales CF24 3AA, UK 

10 Department of Physics, University of Rome "La Sapienza", Piazzale 
Aldo Moro 5, 1-00185 Rome, Italy 

1 1 Department of Physics and Astronomy, University of Pittsburgh, 
Pittsburgh, PA 15260, USA 

12 Code 553/665, NASA/Goddard Space Flight Center, Greenbelt, MD 
20771, USA 

13 Berkeley Center for Cosmological Physics, LBL and Department of 
Physics, University of California, Berkeley, CA 94720, USA 

14 Department of Physics and Astronomy, University of Pennsylvania, 
209 South 33rd Street, Philadelphia, PA 19104, USA 



15 NIST Quantum Devices Group, 325 Broadway Mailcode 817.03, 
Boulder, CO 80305, USA 

16 Department of Physics, University of Toronto, 60 St. George Street, 
Toronto, ON M5S 1A7, Canada 

17 Max Planck Institut fur Astrophysik, Postfach 1317, D-85741 Garch- 
ing bei Miinchen, Germany 

18 Astrophysics and Cosmology Research Unit, School of Mathematical 
Sciences, University of KwaZulu-Natal, Durban, 4041, South Africa 

19 Centre for Astronomy & Particle Theory, School of Physics & As- 
tronomy, University of Nottingham, Nottingham NG7 2RD, UK 

20 Department of Physics, University of Miami, Coral Gables, FL 
33124, USA 

21 Instituto Nacional de Astrofisica, Optica y Electronica (INAOE), To- 
nantzintla, Puebla, Mexico 

22 Columbia Astrophysics Laboratory, 550 W. 120th St. Mail Code 
5247, New York, NY 10027, USA 

23 Institute for the Physics and Mathematics of the Universe, The Uni- 
versity of Tokyo, Kashiwa, Chiba 277-8568, Japan 

24 Institute of Astronomy & Astrophysics, Academia Sinica, Taipei, Tai- 
wan 

25 Department of Physics, University of California Santa Barbara, CA 
93106, USA 

26 Kavli Institute for Cosmological Physics, Laboratory for Astro- 
physics and Space Research, 5620 South Ellis Ave., Chicago, IL 60637, 
USA 

27 Centre for High Performance Computing, CSIR Campus, 1 5 Lower 
Hope St. Rosebank, Cape Town, South Africa 

28 Department of Physics and Astronomy, Haverford College, Haver- 



2 



ACT Collaboration 



1. INTRODUCTION 

Over the past two decades, precision measurements of the 
Cosmic Microwave Background (CMB) have led to remark- 
able advances in our understanding of cosmology. Combined 
with other observations, they have produced constraints on 



TABLE 1 
Scan parameters 



mode ls of the U niverse to percent level accuracy (e.g., Ko- 
|matsu et aL"1 ( |201 Primary CMB anisotropy has been mea- 
sured to cosmic variance p recision by WMAP up to multi- 
poles of approximately 500 (Larson et al. |2011[). Observations 
at finer angular scale s (Friedm an et al.|2 009; Reichard t et al. 
2009b lal IVeneziani et al 



Season 

Elevation 
Scan Width 
Period 
Speed 
Max. Accel. 
Data file length 



2007 

50? 5 

9?6 

19.4s 

l.Odeg/s 

8.1deg/s 2 

15 min 



2008 

50.° 5 
7°0 
10.2 s 
1.5 deg/s 
3.3deg/s 2 
15 min 



2009-2010 

50?5 
7?0 
10.2 s 
1.5 deg/s 
3.3 deg/s 2 
10 min 



2009; 



Kessler et al. ||20 09; Sievers 



etal 120091 |Sharp et aL I 26lOI|Shirokoff et : al.|26ll[|Das et al. 



201 lbt |Reichardt et al.|20TT) , corresponding to the damping mological parameters in Dunkley efaT] ( |2011) . Clu ster de- 



scale of the anisotropics, have led to even tighter parame 
ter constraints, which are improving with each new data set. 
Additionally, CMB measurements at smaller angular scales 
are sensitive to the contribution of point sources, which help 
to re veal the nature of early galaxies or active galactic nu- 
clei ( |Vieira et aLl |2010; Marr iage et al.||2011a) , to the ther- 
mal Sunyaev-Zel'dovich (SZ) effect, which can p rovide inde- 
pendent constraints on cosmological parameters (ISunyaev &| 

~T27TTTf 

d 



Zel'dovic hiri970l [M arriage et al. 2011b; Sehgal etal 
Planck Collaboration VIII| |201 1[ [Reichardt et al.||201 



and 



to gravitational lensing effects on the CM B (Seljak 199 6||Zal-| 
darria gT&Seljak|1999|pas et al.|201 lat|Sherwin et al.|201 1) . 

The Atacama Cosmology Telescope (ACT) is located at 
22°57'35"S, 67°47'13"W on Cerro Toco at an altitude of 
5200 m in the Atacama Desert of northern Chile. Its main pur- 
poses are to map the millimeter wave sky at arcminute scales, 
sampling multipoles up to I ~ 10 4 and detecting and charac- 
terizing foregrounds, including galaxy clusters through their 
SZ signature, and millimeter galaxies. 

Between 2007 and 2011, ACT was equipped with the 
Millimeter Bolometric Array Camera (MBAC), which ob- 
served simultaneously in three bands: 148 GHz, 218 GHz and 
277 GHz. The bands were chosen to avoid major atmospheric 
emission lines and to sample the SZ decrement, null and in- 
crement. Each band had a dedicated set of optics and a detec- 
tor array co mposed of 1024 pop-up TES bolometers ( |Benford| 
|et al.||20"03"]l coupled to the optical signal, called "live detec- 
tors, plus 32 "dark detectors," which w ere not coupled to th e 
sky. For details about the instrument see|Swetz et al.|(|201 1). 



Since first light on October 23, 2007, the telescope has 
had four observing seasons producing more than 90 TB of 
data. Here we describe the 1 1 TB of data taken at 148 GHz in 
2008, including the techniques used for data characterization, 
as well as the data reduction process used to obtain the final 
maps. The 218 GHz and 277 GHz data will be described in a 
later paper. However, a similar data reduction pipeline is used 
to analyze the full dataset. 

The maps obt ained from these data have an angular resolu- 
tion of 1.37' ( Hincks et al.|2010"| and a noise level that ranges 
between 25 and 50 /xK-ar cminT The ca libration of the data 
to WMAP is discussed in Hajian et al. ( 201 1 1. The power 



spectrum of the ma ps is presented in |Fowler et al. ( 2010[ ) and 
Das et al. (201 lb i, with corresponding constraints on cos- 
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tections through their SZ signatur e are p resented by Hincks 
|et al.| ( [20T0) l and Marriag e et al.|(|201 lb|), wh i le extra galactic 
source detections are given in Marriage et al.| ( |201 la|>. Multi- 



wavelength follo w up for these clusters is presented in Menan- 
teau et al. ( 20 1 0)> as well as the discovery a massive cluster at 



z = 0-87, HMenanteau et"aLl[20T2l El Gordo). Their c osmo- 
logical interpretation is discussed in |Sehgal et al. ( 201 l| l. The 
first direct detection of gravitational lensing of the mic rowav e 
background was made using these maps in Das et al.] ( |201 la[ ). 
This in turn demonstrates for the first time that microwave 
background data on t heir own favor cosm ologies with an ac- 
celerating expansion ( |Sherwin et al.|20fT I. 

This paper is organized as follows. A summary of ob- 
servations is given in Section [2] In Section [3] we provide 
background for understanding the sky signal as recorded in 
the time-ordered data stream (TOD). In Sections [4j [5] and [6] 
respectively, we characterize the atmospheric, detector and 
systematic noise found in the TOD. Data calibration into 
units of sky temperature fluctuations, <5Tqvib > is described in 
Section 7 Section [8] describes the pointing solution, while 
Section 9~ explains tne detector time-constant determination 
method. E)ata selection is described in Section 10 provid- 



ing final statistics on the amount of data used for making the 
maps. The mapmaking method is discussed in Section[TT] We 
conclude in Section 12 and present a step-by-step summary of 
the data pipeline from raw data to maps. Appendix[A|provides 
further details about finding and removing correlated modes 
from the data. 

This paper accompanies the public release of the data 
through the NASA's LAMBDA siteQ 

2. OBSERVATIONS 

ACT observes the sky by scanning the telescope in azimuth 
at a constant elevation of 50? 5 as the sky moves across the 
field of view in time, resulting in a stripe-shaped observa- 
tion area. With this scan strategy, the instrument observes 
through a constant air mass, the cryogenics remain stable, the 
telescope's shape remains constant, and the local environment 
and instrumental offsets are sampled in a consistent way. The 
time constants of the detectors, together with mechanical fac- 
tors, limit the scan speed and turnaround acceleration. Table[T] 
gives a summary of the scan parameters used in all seasons. 
The lower acceleration introduced in 2008 was needed to re- 
duce vibrations in the optics (see more details in §|6}. 

The observations are repeated at complementary central az- 
imuth angles to capture both rising and setting skies. This 
cross-linking technique helps minimize systematic effects due 
to the scan in the mapmaking process and improves the deter- 
mination of CMB modes parallel to the scan direction. 

32 http://lambda.gsfc.nasa.gov/product/act. 
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Fig. 1 . — Example TOD from one detector at 148 GHz from October 
21, 2008. This was a good observing night with a PWV of 0.22 mm. 
The slow drift is dominated by changes in the atmosphere brightness. 
The high-frequency noise is dominated by detector noise. Units are 
mK in CMB equivalent units at 148 GHz. The plot displays 3.6 x 
10 5 samples at an interval of 2.5 x 10~ 3 seconds. The telescope was 
scanning while these data were taken. 



Each detector is sampled at a rate of 398.72 Hz and the 
data are stored in 15-minute long data files. They are then 
merged with the rest of the housekeeping data, which in- 
clude the azimuth and elevation encoder readings and time 
of day. The data are compressed to one-third of their origi- 
nal size for storage, using a lossless compression algorithm 
(SLIM)p^l A sample TOD is shown in Figure [l] for a night 
with gooa weather conditions, meaning that precipitable wa- 
ter vapor (PWV) remained below 1 mm. The PWV is a good 
indicator of 148 GHz opacity and overall data quality. 

Despite some variation over the season due to changes in 
sunrise and sunset times, a typical day of observations at the 
site was as follows. The cryogenic system was recycled ev- 
ery day before providing roughly 14 hours of observing time. 
The cryogenic cycle began around 1 1 :00 in the morning lo- 
cal time and lasted for 9.5 hours, so that MBAC was cold by 
about 20:30. At 20:40, warm-up movements for the motors 
and gears were run, and at 20:50 the detectors were b iased 
and the first detector calibration data were obtained (see §7.1| i. 
The observations started at 21:00, usually by scanning the ris- 
ing sky. Around 2:30 more detector calibration data were col- 
lected and the scan was shifted to the west to observe the same 
region that was previously observed while rising. Observing 
in the west has an additional advantage from a hardware safety 
perspective: in the event of a telescope failure, MBAC would 
be pointing away from the Sun at its rising. Observations nor- 
mally ended around 10:40, nearly 2 hours after sunrise above 
the mountains rj At this point, final detector calibration data 
were taken, before the telescope was sent to its home posi- 
tion and the cycle was restarted. When the observable region 
of the sky contained a planet, it was normally scanned every 
other night for calibration and beam measurements. This was 



For further detail see http : / / slimdata . source forge . net/ 
34 For pointing accuracy reasons, we ended up using only data obtained 
less than one hour after sunrise. 



TABLE 2 

Observation summary for seasons 2008 (at 148 GHz). 



Season Dec (Min, Max) RA (Min, Max) Area [deg 2 ] Hours" 
Southern Stripe 

2008 (-57° 15, -48? 1) (20 h 43 m , 7 h 53 m ) 850 1085 

Equatorial Stripe 

2008 (-2? 12, 2? 34) (10 h 21 m , 14 h 48 m ) 280 175 



a Total hours before data selection. 
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Fig. 2. — Histogram showing the PWV during the 2008 sea- 
son, measured at the zenith of the Atacama Pathfinder Experiment 
(APEX) facility. The median value is 0.49 mm. 



not done every single night to avoid producing a gap in the 
CMB map. All of the operations listed above were automated 
and could be performed remotely. 

The 2008 season began on August 1 1 and ended on Decem- 
ber 24, with a total of 136 available nights and 124 nights with 
successful observations, resulting in 1423 hours and 25.7 TB 
of total data for all three frequencies. From the total observed 
nights, four had bad weather conditions, leaving 120 nights 
with usable science observations. The overall calendar time 
efficiency, including day time, was 44%. The median PWV 
during observations across the season, measured at ze nith of 
the Atacam a Pathfinder Experiment (APEX) facility (Giisten 
|et al.|2006] l, was 0.49 mm. Given that the ACT site is about 
140 m higher than the APEX site, we use a correction factor 
PWVact = 0.88PWV A pex to estimate the PWV at the ACT 
site. Figure [2] shows a histogram of the PWV during the sea- 
son. 

Observations were made in two areas of the sky: The equa- 
torial stripe, centered at a declination of 0° (60° and 300° in 
azimuth at 50°5 elevation), and the southern stripe, centered 
at a declination of -53° (150°and 210°azimuth at 50.°5 el- 
evation), covering a wide range in right ascension. Table [2] 
lists the boundaries of the observation areas and the number 
of hours available in each area before data selection. 

3. SIGNAL PROPERTIES 

In this section we characterize the expected TOD signal 
from point sources and extended sources, for comparison to 
the random and systematic noise described later in Sections [5] 
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3.1. Point Sources 



For our purposes, we will define a point source as any ob- 
ject with an angular size comparable to or smaller than the 
beam size of the telescope. The beam full-width at half max- 
imum (FWHM), 6 l/2 , is V.37 for 148 GHz. Planets are the 
most important point sources as their high surface brightness 
makes them useful for both calibration and beam measure- 
ments. Distant galaxies can also be approximated as point 
sources and are helpful for the pointing solution. 

As the telescope scans and the sky rotates, a single detector 
traces a zigzag on the sky, sampling the point source every 
time it intersects with the beam. The number of times a point 
source appears depends on the scan period, the beam size and 
the central scan azimuth. A point source appears as a succes- 
sion of "blips" in the time stream. The shape of these "blips" 
is a slice through the telescope point spread function and de- 
pends upon the angular separation between the center of the 
beam and the location of the source. The angular speed of 
the scan on the sky is given by 9 = v scan cos(50?5). The 2008 
scan speed of v scan = l.°5/s implies an angular sky speed of 
57.'2/s. Assuming a Gaussian-like beam of equivalent width 
and neglecting the transit speed of the source, the 3 dB cutoff 
frequency is 

J-3dB - - n ■ UJ 

7r01/2 

This means that the contribution from point sources to the 
TOD is limited to frequencies below 18 Hz, as can be seen 
in Figure [3] 

Given our scan speed, nominal elevation, sky rotation and 
sampling rate, we expect nearly 10 samples per beam on the 
sky at 148 GHz. 

3.2. Extended Sources 

Using CMB simulations we can estimate the expected TOD 
response to extended sky features. Averaging the power spec- 
tra from many such synthetic TODs, we can estimate the 
CMB power spectrum in TOD space. Figure 3 shows the av- 
erage power spectrum of a simulated 148 GHz observation of 
the CMB sky and the average data power spectrum from one 
15-minute file. Because of Silk damping, the CMB has a rel- 
atively sharp cutoff in its characteristic size: unlike the point 
source observations shown in the figure, the CMB has little 
signal at frequencies higher than 10 Hz. 

Given ACT's angular scan speed, the TOD frequency asso- 
ciated with a sky feature of angular size Of is 



TOD - 



20 



f 



(2) 



where the factor of 2 in the denominator comes from consid- 
ering that the angular scale is the size of a positive or nega- 
tive temperature bump, which is half of a wavelength on the 
sky. Given that multipole moment I relates to angular scale 
as I ~ir/9f for scales small compared to the full sky, an ap- 
proximate conversion between multipole and TOD frequency 
is 

27t/tod 27t/ TO d 



lC os(50?5)' 



(3) 



Note that this relation is only a rule of thumb; the actual 
mapping of TOD frequency into multipole space depends on 



the specific scan strategy. This relation shows that the CMB 
power spectrum in the TOD can be shifted in frequency by 
changing the scan rate. As a reference, the TOD frequency 
corresponding to £ = 3000 was 7.9 Hz in the 2008 season. 
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Fig. 3. — Power spectra of signals in the 2008 season compared to 
the data. The solid curve is the average power spectrum from 737 
live detector TODs from one 15-minute 148 GHz data file, during 
which the PWV was 0.5 mm. The rise at low frequencies is the noise 
contribution from the atmosphere. The dashed curve shows the sim- 
ulated CMB signal. The oscillations in the simulated CMB signal 
are due to enhanced power at harmonics of the scan frequency. The 
dot-dashed curve estimates the point source contribution (consider- 
ing only one hit at the beam center) assuming a Gaussian beam. For 
comparison, the thin and thick dotted curves show the expected re- 
sponse for 218 GHz and 277 GHz, with beam sizes of 1 !00 and 0'91 
respectively. The amplitude of the point source power spectra corre- 
sponds to approximately the signal of Saturn. 



4. ATMOSPHERIC EMISSION 

The atmosphere is a strong emitter and absorber at the 
bands of interest, chiefly due to the excitation of the vibra- 
tional and rotational modes of water vapor. For this reason, 
the PWV is strongly correlated to the level of optical load- 
ing on the detectors. The Atmosp heric Tra nsmission at Mi- 
crowaves (ATM) model (Par do et al.|2001 1 provides an esti- 
mate of the loading as a function of the PWV level. At the 
median value of 0.49 mm during the 2008 observations, the 
loading is approximately 0.5 pW, which corresponds to an 
equivalent Rayleigh- Jeans temperature of 6.4 K at the nom- 
inal elevation. 

During the season, the median atmospheric temperature 
drift over 15-minute observations was 0.22 K (Rayleigh Jeans 
equivalent units as measured by MBAC), with lower and up- 
per quartiles at 0.10K and 0.43 K respectively. 

Turbulence induces a spatial structure in the atmospheric 
signal. Acc ording to the Kolmogorov model of turbulence 
(Tatar ski"i|1961| l, the power spectrum of the fluctuations in a 
large 3-dimensional volume is proportional to q r ~ 11 / 3 , where 
q is the wavenumber. The projected signal observed on the 
sky can have either a q~ n ^ dependence, if the wavelength 
is small enough that the turbulence can be treated as three 
dimensional, or a q~ s > 3 dependence, if the wavelength is large 
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compared to the thick ness of the atmospheric layer supporting 
the turbulent motions ( |Church|1995||Lay & Halverson|2000| l. 

Figure [4] shows examples of the atmospheric signature in 
TOD power spectra. Figure [4a] displays the average TOD 
power spectra from four groups of observations with the tele- 
scope not scanning ("stare" observations). When scanning, 
the knee -where the power law meets the white noise level- 
increases in frequency by around 1 Hz and harmonics of the 
scan frequency leak into the spectrum. Before averaging the 
power spectra from different observations (data files), the av- 
erage power spectrum from the dark detectors was subtracted 
to isolate the atmospheric signal from instrumental 1 // noise. 
TODs were binned as a function of the PWV level measured 
by APEX during the period of data acquisition. It is clear 
how the atmosphere signal grows with PWV, shifting the knee 
towards higher frequencies. On the other hand, the power 
law index stays rather constant near the 2-dimensional regime 
value. The average power spectrum from the dark detectors is 
shown for comparison. It is dominated by the thermal fluc- 
tuations of the cryostat. The fact that the "dark-detector" 
power spectrum appears higher than the others in the fre- 
quency range between the knee and lOHzis the result of the 
power subtraction mentioned above. This indicates that the 
1// plus readout noise dominates over the atmospheric plus 
detector noise in that frequency range. The knee frequency 
ranges from 1 to 5 Hz for stare observations depending on the 
weather conditions. 

The departure from a pure power law shown in Fig. [^im- 
plies that the power law index varies with frequency. By fit- 
ting power laws in small frequency ranges, shown as dashed 
lines in Fig. 4a we were able to measure this dependency and 



group the observations as a function of their power spectrum 
slopes at frequencies near the knee. We found that, for similar 
PWV conditions, the slopes can vary significantly, with in- 
dices ranging between the 2-dimensional and 3-dimensional 
regimes, as shown in Figure 4b Low frequencies, which are 



related to large scales in the sky, are dominated by the 2- 
dimensional regime, with a power law index of approximately 
a = -2.4. At frequencies around 1 Hz the power spectrum 
from some observations (blue dots in Fig. [4b| follow a steeper 
power law, suggesting that, under certain weather conditions, 
the 3-dimensional regime dominates at these smaller scales. 
The slopes tend to increase with wind speed. Higher wind 
speeds are expected to shift larger features of the turbulent 
pattern towards higher TOD frequencies. As large scales are 
dominated by the 2-dimensional regime, this would cause the 
opposite effect of reducing the slopes at higher frequencies, in 
contradiction with what is observed. This result suggests that 
higher winds might be associated with intrinsic properties of 
the turbulent layer, such as its width or height, p roducing a 
steeper power law (Chu rch|1995 Lay & Halverson 2000). 

Atmospheric structures larger than 24' -the field of view 
of a detector array- appear as a common mode among de- 
tectors, while smaller features can produce sub-array "corre- 
lated modes," as described in AppendixlA] Given that the dis- 
tance to the turbulent layer is commonlyTess than a kilometer 
(Robs orTet al.|2002) , it becomes out of focus, smearing out 
to scales of nearly 10', in agreement with our optical simu- 
lations. This corresponds to roughly a third of the array size. 
Then, given the scan speed, such signals appear at frequencies 
between 2 or 3 Hz in the TOD. In general, the common mode 
of the detectors is a good estimator of the atmospheric sig- 
nal, but the estimate can in principle be improved by dividing 
the array to account for sub-array atmosphere structure. Our 
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Fig. 4. — Atmospheric signature in TOD power spectrum, (a) Av- 
erage power spectra for four groups of non-scanning TODs selected 
by PWV level. The average power spectrum from the dark detectors 
was subtracted from each group power spectrum and is shown for 
comparison. The dashed lines show a power law fit to each spec- 
trum. The legend indicates the mean PWV and the power-law index 
from the fit for each group. Logarithmic binning was used to reduce 
the noise in the plot, (b) Power law index as a function of frequency 
for the average power spectra from three groups (blue (•), green (T) 
and red (A)) of observations selected by their power law index at 
frequencies close to the knee, all belonging to the third PWV se- 
lected group in (a). Each value was obtained by fitting a power law 
to the corresponding power spectrum in a small frequency range (in 
logarithmic space); the listed frequencies are the mean values of the 
associated frequency range. The error bars show the dispersion of 
the indexes from the members of each group. The same is shown in 
brown (■) for the dark detector average power spectrum for compar- 
ison. 
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TABLE 3 

Median value of the typical sensitivity per detector and 
the total sensitivity for the array at two frequency 

ranges after removing the main 28 modes of correlated 
noise. the units are equivalent temperature for cmb 
fluctuations. 



Range NET lyp [flKy/s] NET tol [fiK^] 

5-20Hz a 786±28 31.7±2.7 

100-120 Hz 1052 ±21 42.7 ±3.2 



TABLE 4 
Noise contribution summary. 



Noise source 


NET 2 [mK 2 s] 


In-band detector noise 


0.30±0.10 


Aliased noise 


0.17±0.08 


Optical loading 


0.15 


Total (typical) noise 


0.62 ±0.04 



a The signal band is < 30Hz so these entries best estimate the instrument 
sensitivity. 



attempts to detect coherent motions of atmospheric features 
across the array, as would appear in a moving frozen-sheet 
model of the turbulent layer, were not successful. Instead, we 
suppress the atmospheric noise by solving for the strongest 
atmospheric modes in the time stream during the mapmaking 
process, as discussed in j fTT| 

5. SYSTEM SENSITIVITY: UNCORRELATED NOISE 

Above the atmosphere knee frequency, the TOD is domi- 
nated by broadband random noise. This noise is generated at 
the detectors and readout circuitry and it is essentially uncor- 
rected among detectors. After de-projecting the correlated 
modes from the atmosphere and systematics (see Section [6] 
and Appendix [A|, this noise can be measured by averaging 
the TOD power spectral density (PSD) over the desired range 
of frequencies. 

The total sensitivity of the array, expressed as the noise 
equivalent temperature (NET), is given by 



NET fnt = 



'JVdM 

E 

. i=0 



1 



NET? 



-1/2 



(4) 



where NET, is the NET of each working detector. Thus, the 
typical sensitivity per detector can be defined as NET typ = 
NET tot yWdet. where A/det is the number of the "effective" de- 
tectors, as defined in j fT0| 

Table[3]lists the typical values of the total and average NET 
for a mid-frequency range (5-20 Hz) and a high-frequency 
range (100- 120 Hz), in equivalent temperature units for CMB 
fluctuations. Uncertainties show the dispersion among 15- 
minute TODs. The noise increase at higher frequencies is 
d riven by intrinsic properties of t he bolomete rs, as described 
in |Marriage| ( |2006| l, |Zhao et aL| ( |2008) > and |Niemack et al.| 
( 2008 ). However, our signal band is limited to below 30 Hz by 
the telescope scan speed and the beam size (see Figure |3J. 
Thus, the mid-frequency range is the best estimate of the in- 
strument sensitivity. 

As given in Table |4] the typical noise variance per detector 
is 0.62±0.04mK 2 s. The total noise can primarily be sepa- 
rated into in-band detector noise, aliased detector and readout 
noise, and photon noise, all of which add in quadrature. 

To reduce the effects of noise aliasing, the detectors are first 
sampled at 15.15 kHz, then a digital four-pole Butterworth 
low-pass filter with a cutoff frequency of 122 Hz is applied, 
and finally the data are resampled at 398.72 Hz. The detector 
noise bandwidth is limited to 8 kHz by a 700 nH inductor in 
series with each TES. 

To determine the effects of optical loading and aliasing, we 
performed a dark test in which we opened the cryostat, put a 



4K (reflective) cover over the detectors, and collected data at 
a variety of sampling rates. The typical noise level in the dark 
was 0.46±0.14mK 2 s, where the uncertainty shows the dis- 
persion among detectors. The units are equivalent to the ones 
in Table|3] After fitting the noise as a function of sampling fre- 
quency, the in-band detector noise yielded 0.30±0.10mK 2 s 
and the total aliased noise contribution was 0. 17 ± 0.08 mK 2 s. 
The latter includes aliasing from both detector and readout 
noise. 

The readout noise is dominated by SQUID noise and 
preamplifier noise. Fully sampled, the readout noise is esti- 
mated to be around 1.2 x lO^mK^s (based on 50 MHz mea- 
surements), which is expected to increase by roughly a factor 
400 when sampled at 15 kHz, reaching «0.05mK 2 s. The 
SQUID noise aliasing was significantly reduced in season 
2010 by reducing the readout bandwidth. Taking all this 
into consideration, we estimate that slightly more than half of 
the total aliased noise contribution is aliased detector noise, 
which is consistent with the detector noise aliasing analysis 



presented in Niemack (2008]). The other half is SQUID noise. 

The dark tests also revealed that the optical loading con- 
tributes 0.15mK 2 s of photon noise. By measuring the satu- 
ration power of the detectors at different atmospheric condi- 
tions throughout the season, and comparing them to the satu- 
ration power in the dark, we found that the optical loading is 
2.24 pW when the PWV is mm. This loading is dominated 
by spillover contributions, emission from optics and dry atmo- 
sphere emission. The spillover contribution was reduced by 
0.36 pW in season 2010 by adding a baffling structure around 
the secondary mirror. In the same way, we determined that 
water in the atmosphere contributes another 0.7pW/mm of 
loading, so in nominal conditions (PWV = 0.49 mm) the total 
optical loading is 2.59 pW. This contributes to the noise as 
photon noise: for a fully incoherent detector coupling to both 
polarizations, the noise in units of power squared per unit fre- 
quency is given by 



NEP « 2hvP+ ■ 



(5) 



where h is Planck's constant, v = 149.2 ±3.5 GHz is the 
central frequenc}p^j Ai; = 18.4 GHz is the band width, and 
P = 2.59 p W is the power absorbed by the detector ( |Zmuidzi^] 
|nas|2003) . The first term in the equation corresponds to pho- 
ton shot noise while the second term corresponds to noise 
from photons arriving in bunches, as is expected for thermal 
radiation when the occupation number is high. This leads to 
a photon NEP of 3.0 x 10- 17 W/%/Hz, or roughly 0.17mK 2 s 
in CMB temperature units, in agreement with our measure- 
ments. 

35 The central frequency for a Rayleigh-Jeans source. 
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Fig. 5. — Total array NET during the 2008 season in the mid- 
frequency range (5-20 Hz) for 148 GHz. The values are grouped 
in days and the error bars are the standard deviations within each 
day. The PWV is also plotted for reference (dashed line). Values 
are calibrated to CMB equivalent units. Before computing the sensi- 
tivity, eight multi-common modes, four row-correlated modes, four 
column-correlated modes, and the residual twelve modes with high- 
est singular values were removed (see Appendix |A| for details). The 
noise improvement after September 24 came from turning off some 
oscillating detectors, which were contaminating neighbors. 



Table [4] shows a summary of the different uncorrelated 
noise contributions that determine the system sensitivity. 

Figure [5] shows the mean noise between 5 and 20 Hz for all 
data files in the 2008 season. This is the total noise which 
goes into the maps and can be reduced only by increasing the 
o bservation time. Noise e stim ates from the ma ps can be found 
Marriage etaLl ( |2011a| l and |Das et aT1 ( |2011b| >. 



m 



6. SYSTEMATIC NOISE 

In addition to the atmospheric and random noises, the ce- 
lestial signal is contaminated by several systematic noise 
sources. The important ones are thermal drift of the cryo- 
stat, mechanical accelerations which couple both optically 
and thermally (by causing detector temperature oscillations), 
electromagnetic pickup and magnetic pickup. All of these 
effects cause zero-lag correlations among different detector 
TODs, or can be directly measured using the dark detec- 
tors. Detectors are distributed in 32 "columns" and 33 "rows," 
where the last row contains the dark detectors. Each column 
shares the same time-multipexed SQUID readout circuit ( fde| 
|Korte et "aI7||2003) > within which each row is read simultane- 
ously. Thus, there is one dark detector per readout circuit, 
serving to assess systematic noise from it. Moreover, the high 
redundancy of live detectors across the array can also be used 
to assess systematic effects (see Appendix |A). 

6.1. Thermal Drift 

The current signals from the bolometers are amplified in 
100-SQUID series arrays (SA) operated at 3K ( |Swetz et aE] 
201 1) . Slow temperature changes in the SAs produce slow 
drifts in the TODs. During the first 10 hours of the night, 
the SA temperature drifts down by 250 mK, and rises up by 
150 mK in the last 2 hours. In terms of equivalent sky temper- 



ature at 148 GHz, this corresponds to a drift of nearly 4 Kin 
signal. In frequency space, the drift imprints a 1// signature 
on the data, which meets the detector noise level at a knee 
of nearly 1 Hz. Despite small differences in the responses of 
different SA amplifiers, most of this signal appears as a com- 
mon mode to all the detectors. As the coupling occurs at the 
readout circuit, this signal is also present in the dark detectors, 
simplifying its identification and eventual removal. 

In most cases this signal is subdominant to the atmosphere 
signal, but when the PWV is low enough the thermal drift 
becomes comparably significant. This is not evident in Figure 
4a since the examples of low-atmospheric power get averaged 
clown when grouping many TODs. 

In contrast, drifts in detector temperature cause no measur- 
able effect because it is servo-controlled to within less than 
1 mK error. 

6.2. Mechanical Accelerations 

Mechanical accelerations, occurring mainly at the scan 
turnaround, can couple optically and thermally to the detector 
and result in an undesired instrumental response. The opti- 
cal coupling can be produced by a mechanical motion of the 
detector coupling layer, which is a 40 /im thick silicon layer 
placed 100 /im away from the detectors for optical impedance 
matching. The coupling efficiency is sensitive to the distance 
between the coupling layer and the detectors, so small vibra- 
tions can cause significant effects, especially at higher fre- 
quency bands where the coupling layer is thinner and the cou- 
pling efficiency is more sensitive to changes in the distance. 
The effect is also larger near the center of the detector array, 
presumably because the coupling layer vibrates in its funda- 
mental mode. In the TOD, vibrations manifest themselves 
as a series of spikes visible in the common mode at the scan 
turnarounds, with opposite signs for opposite directions, lead- 
ing to lines in the power spectra at several harmonics of the 
scan frequency. In the 148 GHz "waterfall plot" of Figure 
|6j the vibrations contribute to the scan harmonics and some 
resonant lines at higher frequencies. The latter are also seen 
for stare observations, and their central frequencies are shared 
among all three arrays, suggesting that they correspond to nat- 
ural frequencies of the whole system. To mitigate the mechan- 
ical effect, the turnaround acceleration was reduced from its 
value in the 2007 season to the value shown in Table Q] af- 
ter October 8, 2008. The coupling layers for 218 GHz and 
277 GHz were removed after the 2008 season. 

Thermal perturbations of the cryostat, revealed in spectral 
analysis of the bath temperature, also add to the low frequency 
scan harmonics seen in Figure [6] This effect is column- 
dependent and differs from the coupling layer effect in that 
the signal is not stronger at the center of the array. 

In general, high frequency spectral lines in the TOD from 
different detectors destructively interfere when projected into 
map space, mostly canceling out when many observations are 
combined. On the other hand, low frequency harmonics of the 
scan may produce non-negligible bar-like features in the maps 
perpendicular to the scan direction, if not treated properly. 

6.3. Electromagnetic Pickup 

Electromagnetic pickup couples to the readout circuit, pro- 
ducing various signatures in the data. These signals ap- 
pear correlated among subgroups of detector TODs, par- 
ticularly among detectors from the same column or row. 
We call these "column" or "row" correlations respectively, 
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(a) Before removing correlated modes. (b) After removing correlated modes. 

Fig. 6. — Frequency-space waterfall plot containing the power spectrum from all live detectors from the 148 GHz band, before and after 
removing 20 correlated modes. This TOD was obtained under nominal scanning conditions (f sca n = 0.1 Hz), on November 20, 2008. The data 
are calibrated to CMB temperature equivalent units. Each horizontal line in the 2D plot represents the power spectrum of a single detector. The 
black horizontal lines separate detectors from different columns. Most of the low frequency lines are scan harmonics. They are explained by a 
combination of effects including spatial variations in the atmosphere brightness through the nearly triangular scan pattern, thermal oscillations 
of the cryostat and coupling layer oscillations. The thick lines around 5 Hz (see columns 2-5) are also present in stare observations. They are 
thought to be resonant frequencies of the system. Note that the science band is between about 1 and 20 Hz. After removing 20 correlated modes 
from every detector TOD, harmonic features are significantly reduced, as is the low frequency power from the atmosphere and from thermal 
drifts. When making maps, the correlated modes are identified and properly de-weighted to produce an unbiased map solution. 



each corresponding to a different source of electromagnetic 
pickup. Narrow-band signals appear correlated among detec- 
tors within the same column. They couple in somewhere af- 
ter the first SQUID stage of the readout circuit, as the sub- 
sequent circuitry is shared by the detectors in a column. On 
the other hand, broad-band signals appear correlated among 
detectors within the same row, or even a few rows apart. 
We believe these are caused by rapidly varying signals, such 
as spikes from strong current transients: given our time- 
multiplexed readout scheme, in which the same row is read 
from all columns simultaneously, these signals become corre- 
lated among detectors in the same row. 

Figure [7] shows the detector correlation matrix for a given 
15-minuteTOD file, where every element in the matrix corre- 
sponds to the correlation between the TODs of detectors i and 
j, which index the elements of each axis. For instance, the 
diagonal elements correspond to i = j, so they are all equal 
to one. In Figure [7a] adjacent detectors on each axis belong 
to the same column in the array, with a black line separat- 
ing detectors from different columns, while Figure 7b is or- 
ganized by rows. We can clearly see the kinds of correlations 
described above affecting either columns or rows of detectors. 
As would be expected from a stochastic distribution of tran- 
sients, the patterns seen in the row-correlated matrix change 



in time. To quantify the correlations, a quality factor is de- 
fined as the mean of the squared off-diagonal elements of the 
correlation matrix: 



Q- 



N(N-l) 



N 



(6) 



where Cij is the correlation between the TODs of detectors i 
and j, and N is the number of detector TODs. Lower Q factors 
indicate less correlated noise. 

The column-correlated signals were significantly reduced 
by carefully shielding the system from electromagnetic noise 
in the environment. The row-correlated signals, instead, were 
found to be related to the switching power supplies feeding 
the readout electronics. They were replaced by linear power 
supplies at the end of the 2008 season. 

If not treated properly, broadband row-correlated signals 
can produce features in the maps, mostly seen as features per- 
pendicular to the scan direction for that particular observation. 

6.4. Magnetic Pickup 

The readout circuit of the detectors uses SQUIDs, which 
are very sensitive to magnetic fields. Despite significant mag- 
netic shielding, some residual pickup remains as the telescope 
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(a) Column dominant. Magnetic pickup and atmosphere removed. 
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(b) Row dominant. Magnetic pickup and atmosphere removed. 
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(c) Column dominant. Magnetic pickup, atmosphere and dark modes removed, (d) Row dominant. Magnetic pickup, atmosphere and dark modes removed. 

Fig. 7. — Correlation matrices for a 15-minute data stream obtained December 7, 2007. The 551 active detectors are each numbered, starting 
from the upper left of the correlation matrix, in order of columns (left panel s) or rows (right panels). The remaining 473 detectors were cut 
from this particular data stream due to inadequate performance (see Sec. 1 1 0. 3| > - The solid lines divide the 32 individual columns (left panels) or 
rows (right panels); the different widths between the solid lines reflect the varying number of active detectors in each column or row (cf. Fig. [8}. 
The magnetic pickup and the atmospheric signal have been removed from the time stream by subtracting a best-fit sinusoid for the former and 
by removing eight multi-common modes for the latter (see Appendix[A](. The top panels show these filters only, with a quality factor Q = 0.081 
(Eq. |6j. The bottom panels additionally have dark correlated modes removed, with an obvious large suppression of detector correlation and a 
quality factor of Q = 0.015. 



sweeps through the Earth's magnetic field. Although the scan 
is almost triangular, the observed signal is chiefly sinusoidal, 
which can be explained by eddy current losses and hystere- 
sis in the magnetic shielding. The magnetic signal amplitude 
and phase are fairly stable across detectors in the same col- 
umn, but differ significantly between columns. We believe 
these phase shifts are related to the complexity of SQUIDs 
and electrical loops in the system. A typical amplitude of this 
signal is equivalent to nearly T ~ 7mK in CMB units. 
As this is a readout-related signal, it can be measured using 



the dark detectors and removed by fitting and subtracting a 
sinusoid. Moreover, the scan frequency is significantly below 
the science band and the magnetic pickup is heavily down- 
weighted in the mapping. 

7. CALIBRATION 

Our calibration takes into account the detector properties 
and electronics; camera and telescope optics; and atmospheric 
conditions. It consists of three steps. The first adjusts for vari- 
ations in detector responsivity due to changing atmospheric 
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loading and atmospheric opacity. Next, a detector-detector 
gain term is applied to account for relative variations in opti- 
cal coupling across the camera. Finally, a single normaliza- 
tion is applied to all the data to account for the overall system 
efficiency for celest ial sources. Th is procedure is similar to 
the one presented in Switzer (2008). 

7.1. Calibration Variations with Time 

The nightly variations in atmospheric loading and daily 
cryogenic cycling necessitate a rebiasing of the detectors at 
the beginning of each night's observations, which affects the 
detector response. Additionally, the opacity of the atmosphere 
varies from night to night. Changes in loading and opacity 
throughout the night produce smaller variations in the system 
response. 

In order to account for the variations in system response 
between nights, a responsivity calibration is derived from the 
analysis of /- V curves taken during the nightly rebiasing. An 
I-V curve is the relation between the output current response 
(7) of the detector's SQUID readout circuitry to a slowly 
ramping detector bias voltage (V). For the 148 GHz band, the 
median deviation in responsivity between nights as measured 
by the I-V curves is 3.0%. 

A second test, known as a bias step ( Niemack]200 8), is per- 
formed three times per night to detect changes in detector re- 
sponsivity through the night. This is achieved by recording 
the detector response to a series of small, square-wave pulses 
applied to the detector bias voltage. The responsivities can 
be obtained by analyzing the detector responses to this signal, 
and are found to be in agreement with responsivities obtained 
from the I-V curve analysis. For the 148 GHz band, the me- 
dian deviation in responsivity over a night, as measured by the 
bias step, is 1 .0%. For more information about the ACT detec- 
tors, the biasing routin e and responsivity, se e |Fisher| ( 2009 ), 
|Battistelli et~aL] ( |2008] >, |Swetz et al.| pOTT) and |Zhao et at] 



from |Griffin & Orton| ( [T993] l, in combination with WMAP-7 
measurements of the brightness of Mars and Uranus ( Weiland 



et al.|201 1)>. It yields 103.5±6Kat 148 GHz, producing a 6<% 



(2008). The time dependence on the atmospheric conditions 
will be discussed below together with the overall system cali- 
bration. 

7.2. Time Independent Relative Detector Calibration 

To determine the detectors' relative gain coefficients, we 
use the large common-mode signal provided by the variations 
in atmospheric brightness. This is analogous to a flat-field cal- 
ibration done with an optical CCD. For each 15 -minute data 
file, we compute the gain factor that best fits the detector drift 
to the common mode drift. The dispersion of this factor per 
detector over the season, averaged over all detectors, is better 
than 2% rms for the 148 GHz array. This includes the variabil- 
ity of the time-dependent calibration step. We correct for this 
by multiplying each detector time stream by its corresponding 
gain factor. 

7.3. Overall System Calibration 

An overall system calibration is established by measure- 
ments of Uranus. By comparing the peak response of the 
instrument to the known brightness temperature of Uranus 
(taking into account the ratio of Uranus' solid angle, f2u = 
0.235 ±0.010nsr, to the telescope's solid angle at 148 GHz, 
^act = 215.8 ± 1.5nsr), an overall calibration of the system 
response is obtained. The statistical error derived from 30 
measurements in the 148 GHz band is 1%, while solid angle 
uncertainties are below 1%. The equivalent Rayleigh- Jeans 
temperature of Uranus is obtained by reprocessing the data 



uncertainty in the overall calibration. 

The conversion to sky temperature also depends on the at- 
mosph eric transmission. This is estimated with the ATM 
model ( |Pardo et al.||200l) , which uses PWV measurements 
from APEX corrected for the ACT site and other Atacama- 
specific parameters. When PWV measurements are not avail- 
able, the season-average transmission (T = 0.976) is used. 
For the 148 GHz band, the rms of the transmission during the 
2008 season was 2.3%. Given that large atmosphere loading 
can excite non-linearities in the detector responsivities, an ex- 
tra degree of freedom is given to the planet fit as a function of 
PWV, producing a final correction to the calibration factor. 

The overall map calibration is compared to the WMAP-7 
year map, by correlating multipoles in the range 400 < I < 



1000, providing an uncertainty of 2% in temperature (Hajian 



|et al.|20ri) . This analysis shows that the Uranus calibration 
is a factor of 6% lower than the WMAP calibration, which is 
within the expected uncertainty. 

Putting all pieces together, the overall system calibration 
for 148 GHz in the 2008 season is 

F(w) = C exp [{T d + r w w+t x (w- w)) I sin 6> a i t ] , (7) 

where C = 19.41 K/pW is the overall calibration factor, 9. d \ t 
is the observation altitude, Td = 0.0093 is the "dry opacity," 
t w = 0.0190mm" 1 is the "wet opacity," r A = 0.0138mm" 1 is 
the the extra degree of freedom, w is the PWV measurement 
from APEX, corrected for the ACT site in millimeters and 
w = 0.44mm is a pivot PWV used in the fit. Here tj and r„. 
are fixed parameters obtained from the ATM model, while C 
and t x are free parameters. If t w is allowed to vary, we recover 
the sa me value that is predicted by the ATM model (Hincks 
[2009l l. 

8. POINTING 

The pointing solution is decomposed into relative pointing 
between detectors and absolute boresight pointing, both of 
which we discuss here. 

8.1. Relative Detector Pointing 

The relative pointing between detectors was determined by 
modeling the beam in the time streams of planet observations. 
This analysis made use of approximately 30 observations of 
Saturn at the normal CMB observing altitude of 50?5, in two 
azimuth ranges corresponding to the rising and setting of the 
planet. 

The telescope scans are rapid enough that each detector 
samples the vicinity of the peak response to a planet several 
times in a single observation. It is thus possible to model the 
beam position and shape in two spatial dimensions, project 
this to a detector signal using the telescope encoder informa- 
tion, and fit the data in the time domain. For simplicity, we 
use a 2-dimensional Gaussian as the beam model. The fit pro- 
duces azimuth and altitude offsets for each detector relative 
to the telescope pointing encoders, as well as measures of the 
peak response, beam FWHM, and optical time constants. 

When combining different planet observations, we first 
aligned them by applying an offset correction to each one. 
The position of each detector is taken as the mean of the posi- 
tions from all observations, after rejecting outliers. This pro- 
duces relative detector positions for each array that are in good 
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agreement with design expectations ( |Fowler et al.||2007 i. A 
comparison between the offsets for rising and setting observa- 
tions (which differ in azimuth by approximately 95°) shows 
no significant rotation or shear of the array relative to the lo- 
cal altitude and azimuth axes. This constrains the tilt in the 
telescope azimuth axis to be smaller than 1' and indicates that 
rotation of the array with respect to azimuth is negligible. 

The uncertainty in the relative pointing is no greater than 
1"5 for all three arrays. Given the large number of detectors 
in each array, this small uncertainty in the relative detector 
positions amounts to a negligible contribution to the pointing 
error in the final maps. 

8.2. Absolute Boresight Pointing 

The azimuth and altitude of the telescope encoder readings 
must be corrected to account for their offset from the true 
boresight for each frequency band. The correction is differ- 
ent for each of the four central CMB observation azimuths, 
namely the equatorial and southern stripes at rising and set- 
ting orientations. 

The correction is obtained by projecting the data from a par- 
ticular band and orientation into a map (after having applied 
the relative pointing solution described above), and compar- 
ing the positions of the bright point sources to catalogue posi- 
tions. The resulting offsets in equatorial coordinates are con- 
verted to offsets in the boresight position, and this procedure 
is iterated to ensure convergence. The boresight offset varies 
by about l'over the season, primarily in elevation. 

The dominant source of error in the offset correction comes 
from estimating the centroid positions of the point sources in 
the maps, before matching them to the catalog positions. This 
error scales inversely with the square root of the number of 
point sources available, which was 20 for the 148 GHz band, 
and with their signal to noise ratio. The uncertainty was found 
to be 2 "6 for the southern stripe and 5 "3 for the equatorial 
stripe at 148 GHz. These values were obtained by adding in 
quadrature the error in the fit from the rising and setting maps. 

The telescope pointing is expected to vary slightly due to 
thermal deformation of the mirror structure. This variation 
is estimated from observations of Saturn taken at nearly the 
same azimuth and altitude on different nights. The rms point- 
ing variation of such observations is 4."3. Since any pixel in 
the final season maps contains contributions from many dif- 
ferent nights, this random variation does not contribute signif- 
icantly to the pointing uncertainty. 

Pointing deviations are significantly higher at dawn when 
the telescope temperature begins to change more rapidly, 
showing a trend that repeats every morning. The drift begins 
nearly 1 hour after sunrise, reaching nearly 50" an hour later. 
Data taken more than 1 hour after sunrise are not included in 
the maps. 

On top of this, the altitude was observed to drift by about 
20" over the course of the 120 day season, producing pointing 
errors that were correlated with right ascension. To remove 
this trend, linear corrections of 0"2/day and 0"15/day were 
applied to the rising and setting fields respectively. These cor- 
rections have little effect on the maps other than to remove the 
small, RA-dependent pointing offset. 

The net pointing uncertainty is thus dominated by the sys- 
tematic uncertainty in the alignment of the rising and setting 
maps, along with some residual variation due to temperature- 
dependent mirror deformation. Comparing the positions of 
bright point sources in the maps to their catalog positions, 
we estimate the pointing error in the final maps to be 4"8, 



with no preferred orientation (Marriag e et al.| 2011ai. Note 
that this uncertainty is much smaller than the beam size, and 
would cause errors of less than 1% in the measured flux for 
point sources. 

9. DETECTOR TIME CONSTANTS 

A detector's time response is limited by its electro-thermal 
properties. We model their optical step response as an expo- 
nential decay with a time constant r. A finite response time 
results in a small shift in the spatial position of a point source. 
The shift depends on the scan direction. In the fits for the rel- 
ative pointing solution described above, time constants were 
included in the beam model as a single-pole low-pass filter 
in the time domain. The time constants are measured with 
a precision of St < 0.5 ms, which is an upper bound to their 
dispersion from the analysis of 30 Saturn observations. The 
median time constant of the 148GHz array is 1.9 ± 0.2ms 
(hdB = 83.8 ± 8.8 Hz), with only a handful of detectors show- 
ing responses slower than r = 10ms (f^B ~ 15Hz). For com- 
parison, Equation[3]implies that I = 10 4 corresponds to 27 Hz. 

10. DATA SELECTION 

Data selection can be divided into two types: data file selec- 
tion (all detectors) and single-detector TOD cuts. The former 
determines the number of observing hours, while the latter de- 
termines the number of "effective detectors" within each data 
file, defined as the sum of the fraction of the time that each 
detector passed the data selection: 



'total 



(8) 



where 7} uncut is the available time for detector i after data se- 
lection, and 7i t a i is the total time before data selection. The 
number of array-wide effective detectors can then be defined 
as the average of N e s over all the available data files in the 
season. 

The remainder of this section presents the data selection 
methods and statistics for these two types of data cuts. 

10.1. Detector Classification 
The detectors are classified in three groups: 

• Live Detector Candidates: Those that have the potential 
to be used for map making; 

• Dark Detector Candidates: Those that do not couple to 
the sky (for example a broken pixel) but can be used to 
diagnose systematics. This includes the subset of the 
32 dark detectors which work properly and the subset 
of defective detectors with working readout circuit; 

• Broken Detectors: Those defective detectors that can- 
not be used to probe systematic effects. They include 
those former live and dark detectors with defective 
readout circuits, and slow detectors (r < 10 ms). Many 
of them must be turned off while observing to prevent 
them from interfering with other detectors. 

Several methods were used to classify detectors into these 
categories. These include assessing correlation with the atmo- 
spheric signal (by far the la rge st signal), searching for consis- 
tently oscillating detectors^] and finding biasing problems. 

36 We have implemented automatic methods to identify oscillating detec- 
tors on the fly and disconnect them. 
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TABLE 5 

Summary of number of detectors in each of the three 
groups, number of detectors cut by each criterion and 
effective number of detectors for 148 ghz in 2008. 
errors are the standard deviation of the number of 
detectors. 



Number 



Detector classification 

Live Candidates 795 

Dark Candidates 128 

Broken Detectors 133 

Total 1056 



Cuts by criterion applied to live detectors 



Calibration 


3±4 


Drift 


95 ±61 


Correlation 


8±16 


Gain 


2±8 


Mid-F. Noise 


41 ±33 


HFRms 


15 ±12 


HF Skewness 


3±4 


HF Kurtosis 


7±9 


Scan 


2±9 


Glitch 


26 ±23 



Effective detectors 

593 ± 95 



Table [5] gives a summary of the number of detectors in each 
of the three groups for each array. 

Once the live and dark detector candidates have been iden- 
tified, a number of possible pathological behaviors may still 
justify removing part of the data from the reduction pipeline. 

10.2. Data File Selection 

Out of the total number of data files acquired, we rejected 
files for the following reasons and in the following order: 

• Files that correspond to planet observations and other 
calibration or engineering tests; 

• Data taken more than one hour after sunrise, to avoid 
pointing and beam errors caused by telescope deforma- 
tions as it was thermally settling; 

• Files with fewer than 400 effective detectors, as they 
were considered likely to be pathological; 

• Bad weather: PWV greater than 3.0mm (transmission 
below 90%); 

• Poor cryogenic performance: Detector base tempera- 
ture more than 7 mK above the nominal temperature, or 
when it changed more than 1 mK within the 15 -minute 
file; 

• Poor calibration: If the relative gain dispersion of the 
detectors was more than 10%; 

• Data for which the analysis software failed. 

The final amount of data available for analysis is summa- 
rized in Table [6] 



TABLE 6 

Data file selection in 2008 at 148 GHz. 



Type 


Obs. Hours 


% 


V-tllCllU.a.1 


3264 h 


100% 


Total observation 


1423 h 


43.6% 


Total survey 


1260h 


38.6% 


Later than 1 h after sunrise 


-198h 


(6.0%) 


Low effective detectors 


-156h 


(4.8%) 


High PWV conditions 


-41 h 


(1.2%) 


Cryogenic problems 


-28 h 


(0.9%) 


High gain dispersion 


-lOh 


(0.3%) 


Other 


-13h 


(0.4%) 


Uncut South 


772 h 


23.7% 


Uncut Equator 


44 h 


1.3% 


Uncut Total Survey 


816h 


25.0% 



10.3. Detector Cuts 

Detectors are affected by sporadic pathologies. Depending 
on the kind of pathology, it may be necessary to reject a sec- 
tion or the full length of a detector TOD from a given data 
file. 

The main causes for these pathologies are quantum jumps 
of the magnetic flux of a SQUID in the readout circuit (V-<fi 
jumps); excessive detector noise; conducted noise from oscil- 
lating detectors; excessive electromagnetic pickup; and me- 
chanical contamination, which can be optically or thermally 
coupled into the signal. 

The following tests were performed over 15 -minute data 
files to detect these pathologies: 

• Drift test: This probes low-frequency deviations of a 
detector TOD from the atmospheric signal. The data are 
first low -pass filtered above 50mHz and calibrated into 
units of power, as described in ^7] Then the thermal 
drift is removed by de-projecting both the dark detec- 
tor common mode and the housekeeping temperature 
of the detectors. Finally, the atmosphere signal is re- 
moved by de-projecting eight multi-common modes, as 
explained in Appendix [A] The drift error is defined as 
the standard deviation of the residual data. Detectors in 
148 GHz are cut as outliers if their drift error is higher 
than 0.35 fW. 

• Correlation test: The Drift test is complemented by 
finding the correlation between the detector drift and 
array common mode drift for TOD frequencies below 
50mHz. Detectors that correlate less than 98% are ex- 
cluded. 

• Gain test: The Drift test assesses the shapes of the 
TODs, but not their amplitudes. The relative detector 
gain can be quantified by the factor that best fits the at- 
mospheric drift to every single TOD at frequencies be- 
low 50mHz. For this the atmospheric drift is estimated 
as the array common mode after removing the thermal 
contamination. Detectors are cut whenever their gains 
differ from the mean by more than 15%. 

• Mid-frequency noise test: Some pathologies, mainly 
associated with mechanical contamination, become 
prominent at frequencies between 0.3 and 1.0 Hz, be- 
low the "science band." To isolate pathological detec- 
tors, we band-pass filter the data below and above those 
frequencies, de-project the array common mode (also 
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filtered), and obtain the standard deviation of the resid- 
uals, which we call mid-frequency error. Detectors are 
cut whenever their mid-frequency error is more than 2.5 
times the median value for the data file. 

• High-frequency noise tests: At frequencies above 
3 Hz the data start being dominated by detector noise, 
which is chiefly Gaussian. Non-Gaussianities above 
this frequency are mostly caused by electromagnetic 
pickup, conducted electrical noise or opto-mechanical 
perturbations. To isolate them, we high-pass filter the 
data below 5 Hz and then test for non-Gaussianities by 
computing the standard deviation, skewness and kurto- 
sis. This is done within sections of the TOD of one 
scan period, with a characteristic le ngth of 10.2 s. Us 
ing the transformation proposed by D'Agostino & Be 



langer ( 1990), the last two statistics yield a normal dis- 
tribution in the case of Gaussian noise, which is veri- 
fied for most of our data. Outliers are identified by how 
much they deviate from the mean. Only sections of the 
TOD with noise rms lower than 0.9 fW are kept. 

• Scan test: Motion defects and encoder errors are also 
detected and cut, affecting sub-sections of the TOD. 
This is done by examining the azimuth time stream and 
searching for interruptions in the scan pattern. 

• Glitch test: The data are also affected by spike-like 
glitches, for example from cosmic rays. Spikes larger 
than 10 times the noise rms are cut, including a 0.5 s 
buffer on either side. Also, whenever two such cuts are 
separated by less than 5 s they are stitched together into 
a single cut. If more that 5 glitches are found in a single 
detector TOD, then the whole detector TOD is cut. 

• Calibration Bolometer test: During 2008 observations 
a calibration bolometer was used to load the detectors 
with radiation of roughly 400 mK every 24 minutes, 
each event lasting for 1.3 seconds. For these, a window 
of nearly 3 s is excised around the event j^] We included 
this within the "Glitch" cuts in Table 

As a general rule, if more than 20% of the detector TOD 
would be cut, then the full detector TOD is cut instead. 

Dark detectors were selected in an analogous but simplified 
way. In this case, only full detector TOD cuts were performed. 
The cut criteria were the drift error of the dark detectors, the 
gain with respect to the dark common mod^j and the noise 
rms, all given in raw data units. 

10.4. Data Selection Results 

After applying all the selection criteria given above, and 
considering only the data files available for analysis, the av- 
erage number of effective detectors was 593 ± 95 in the 2008 
season at 148 GHz. The error here is the standard deviation 
over different data files. 

Table 15] shows a summary with the number of detectors in 
each of the three detector groups for the 148 GHz array, as 
well as the cut contribution from the eight main selection cri- 
teria: drift, gain, mid-frequency noise, noise rms, noise skew- 
ness, noise kurtosis, and partial cuts with more than 20% of 

37 The calibration bolometer was not used after 2008. 

38 Notice that the common mode of the dark detectors is driven by the 
thermal drift of the cryostat, which is the second largest signal in the live 
detector data. 




Fig. 8. — Percentage of time that detectors were cut across the 
148 GHz array. Each small square represents a single detector. They 
are oriented as the array is projected on the sky. Notice that some 
rows and columns are always cut, which is mainly due to problems 
in the biasing (rows) and readout circuits (columns). 




Fig. 9. — Effective number of detectors in the 148 GHz array in the 
2008 season. Circles denote daily averages and the error bars are the 
standard deviation within that day. The PWV level is shown with 
the dashed line. The increase in the number of effective detectors 
starting in September 24 occurred after some oscillating detectors, 
which had been contaminating other detectors, were turned off. 



the TOD cut. The bottom line is the average number of ef- 
fective detectors in the season. Figure [8] is a diagram of the 
148 GHz array showing the fraction of the time that each de- 
tector is uncut. Figure [5] shows the daily number of the effec- 
tive detectors during the 2008 season, along with the weather 
conditions indicated by the PWV. 

On top of these cuts, a small fraction of the remaining data 
is cut during the mapmaking process, as described in the fol- 
lowing section. 
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11. MAP MAKING 

11.1. Mapping Essentials 

To make maps from ACT data, we solve for the best-fit sky 
given the noise in the data. In particular, we find the sky map 
that minimizes x 2 given a model for the noise, and a model 
for what the data should look like: 



d = Mx + n. 



(9) 



Here x is the model for which we wish to solve, M describes 
how the data depend on the model, and n is the particular re- 
alization of the noise in the ACT data. Traditionally, x is a 
vector whose components are the sky map pixels and M is the 
pointing matrix. In its simplest conceivable form, each data 
point sees a single pixel in the map, so each row of M (corre- 
sponding to a single data point) has a single 1 in the column 
corresponding to the map pixel at which it was pointed. How- 
ever, our model for the data, x, can contain many contributions 
in addition to the sky map: components include atmospheric 
noise, correlated signals in the data, point source fluxes, and 
missing (cut) data samples as discussed below. The mapping 
formalism easily generalizes to cover multiple components as 
long as the data depend on them linearly: 



M=[MiM 2 ---M c ] 



(10) 



(11) 



With an estimate of the noise covariance N = (nil 7 "), which 

is a n samp ie by « samp ie matrix (« samp ie ~ 10 9 ), we wish to find 
the model that maximizes the likelihood function 



£ = exp 



-(d-Mx) r N" 1 (d-Mx) 



(12) 



where d is a vector containing all the data samples. The stan- 
dard linear least-squares solution is 



M r N" 



Mx 



:M r N" 



(13) 



The matrix M r N _I M is too big to be practicably inverted di- 
rectly (we typically have 10 7 map pixels), so we instead itera- 
tively solve the least-squares equation for x using a Precondi- 
tioned Conjugate Gradient (PCG) schem e ( Wright et al.|1996[ 
|Hinshaw et al.|[2003| |Press et al.|2007| >. Preconditioning in- 
volves introducing P, an approximate inverse of M r N~'M, in 
order to speed up convergence of classic Conjugate Gradient, 
and solving the better-conditioned system: 



PM r N _I Mjc: 



PM^'d. 



(14) 



To map ACT data, we use the mapmaking code Ninkasi 
and run on the Scinet General Purpose Cluster (GPC) fLoken 



et al.||2010 >. The mapmaking algorithm for the maps in this 



release is improved from that used in previous papers ( Fowler 

ljjbas 



20 

II) in four 



et al.12 010; Marriag e et al.|201 lalb]|Dunkley et al. 
et al.pOl lbj |Hajian et al.||201 l||Sehgal et al ||20! 
ways: 1) rather than solving for detector-correlated noise (in 
eluding atmospheric noise) explicitly, we put the correlations 
in the noise matrix; 2) we explicitly solve for cut data (time- 
stream gaps); 3) we subtract models for the point sources in 
the ACT maps directly from detector time streams; and 4) we 



re-estimate the noise after creating an initial map to remove 
signal-induced bias in the noise estimation. 

The mapping is done in a cylindrical equal-area projection 
with a standard latitude of S = -53° .5 and pixels of 30" x 30", 
roughly one-third of the beam FWHM. 

11.2. Pre-processing 

In addition to data selection, calibration and pointing, there 
are a few other pre-processing steps that are done to the data 
before solving for the maps. 

We remove the median of each detector time stream for 
each 15-minute period, and subtract a single array-wide slope 
across the period so that the ends of the time streams roughly 
line up. We do this to reduce ringing in Fourier transforms 
and to facilitate searching for correlations among the detec- 
tors. We linearly interpolate across gaps in the data, such as 
those arising from cosmic ray hits (again, to reduce Fourier 
artifacts: the data in the cuts are otherwise not used). We de- 
convolve the time streams by the anti-alias filter (described in 
and the detector time constants (discussed in ^9]). 

Next, the non-optical contamination signals are reduced us- 
ing the dark detectors. This is done as follows: we take the 
time streams from each dark detector and subtract a mean and 
a slope. We then high-pass filter the dark time streams, fil- 
tering out any signal below 5 Hz. Most of the remaining sig- 
nal is in only a few independent modes. We take those cor- 
responding to the seven largest eigenvalues of the resulting 
covariance matrix and do a linear least-squares fit to the live 
detector data, which had been similarly-processed (having re- 
moved the mean and slope, and high-pass filtered). Finally, 
we subtract the reconstructed fit from the data. We can do this 
because the modes subtracted are not correlated to the sky 
signal^ 

In [Fowler et al.| (2010), we downsampled the timestreams 
from 400 Hz to 200 Hz, using a time-domain triangular ker- 
nel of the form [0.25 0.50 0.25]. The ACT signal band (nearly 
30Hz according to Sec. [3l is well away from the downsam- 
pled frequency limit. We find that downsampling does af- 
fect the raw power spectrum at up to several percent for 
£ >^ 5000. This should be mostly corrected for when us- 
ing a beam measured using maps derived from downsampled 
data. However, we also find that the SNR on point sources is 
2% higher for non-downsampled maps. Consequently, we do 
not downsample the time streams in t he re sults presented here. 

Also in contrast to |Fowleret al. (20101, we carry out one 
further step in the time domain: With the 30" pixels some- 
what undersampling the ACT 1 .'37 beam, we find that, to get 
percent-level accuracy in source fluxes, we must deal with 
the bulk of the source flux directly in the time streams. To 
do this, we make first-pass maps in which sources are found 
and their fluxes estimated. Simulations show that these fluxes 
are typically accurate to a few percent, with source flux sys- 
tematically underestimated by approximately 4%. We then 
take these source fluxes and model them directly in the time 
streams using the full ACT beam, and subtract the model from 
the time streams. We do a simple source-only projection into 
a map, and add this to the (mostly) source-free maximum like- 
lihood map. We find that, in simulations, this recovers mean 
source fluxes to 1 % accuracy, at which point it is subdominant 
to the calibration uncertainty. We do no additional filtering in 
the time domain. 
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11.3. Noise Model ing 

The noise structure is modeled in both frequency and time 
domains. We include, as a term in the noise, a time domain 
windowing of the first and last 20 seconds of each time stream 
of the form (1 -cosx) /2, again to reduce Fourier ringing. We 
then search for correlations across the data by examining their 
covariance matrices, split at 4 Hz (£ ~ 1500). We find all 
eigenvectors in the low-frequency covariance matrix with cor- 
responding eigenvalues greater than 3.5 2 times the median 
(the eigenvalues of the data covariance matrix A r A are the 
squares of the corresponding data singular values, so this cor- 
responds to an amplitude of 3.5 times the median in the time 
streams). We then project those eigenvectors out of the high- 
frequency covariance matrix, and again find all eigenvectors 
with eigenvalues larger than 3.5 2 times the median in the re- 
maining high-frequency matrix. We find that the maps are 
not particularly sensitive to the exact threshold values chosen. 
This procedure typically finds 15 to 20 low -frequency eigen- 
vectors and one or two high-frequency ones. The eigenvec- 
tors are response patterns of correlations across the array, and 
usually correspond to things like a common mode, gradients 
across the array, and the row-correlated noise. They provide 
the linear combination of detector time streams needed to pro- 
duce acorrelated mode, represented by the vector \ in Equa- 
tion Al This projection corresponds to removing around 20 



modes out of the approximately 600 available. 

With the shapes of the array correlations in hand, we can 
use them to complete the description of the noise. We solve 
for the correlated modes corresponding to the eigenvectors 
and subtract them from the time streams. We then Fourier 
transform the (de-correlated) detector time streams and the 
correlated modes, and model them using frequency bins. In 
each frequency bin, we find the average power in each detec- 
tor time stream and in each correlated mode. If the detector 
noises are denoted by the diagonal matrix No, the detector 
covariance eigenvectors by the ^detector by n mo de matrix V and 
their corresponding bin-wise noise powers by Ny, then the 
bin-wise Fourier-space noise is simply N/ = No + VNi/V r . 
Here N/ is an «detector by ^detector matrix acting on a single 
frequency bin. This matrix can be quickly inv erted using 
the Sh erman-Woodbury formula (Duncan (1944); see |Hager| 
( 1989) for a review). If we denote the time-domain edge ta- 
penng operator by W and the Fourier transform operator by 
F, then our entire noise inverse becomes 



N" 



-1 = W r F r Nf 1 FW, 



(15) 



where the operator N F l is conformed by all the N^ 1 such that 
it acts on every frequency bin in the Fourier transform. In 
the previous equation, every operator can be represented by 
an n S ampie by n S ampie matrix, in which detector operations like 
N/ are expanded such that all samples within a single de- 
tector are treated in the same way. To minimize our sensi- 
tivity to any low-frequency non-Gaussian component of the 
atmospheric noise, we taper the Fourier weights at frequen- 
cies below 0.5 Hz, with the weight explicitly set to zero below 
0.25 Hz. We note that the noise in the mapmaking equation 
(Eq. [T4| can be interpreted as a set of weights. The map is op- 
timal if the weights are perfect, but the map remains unbiased 
as long as identical weights are used for the left and right sides 
of the equation, and the weights are uncorrected with the sky 
signal. This is the essential part of the method. Time-stream 
filters are by design biased and must be accounted for with 
simulations. On the other hand, in our treatment, modes in the 



map that might cause problems are de-weighted as opposed to 
being filtered, producing an unbiased solution through careful 
consideration of the noise structure of the data. 

One furth er difference remains between these maps and 
those used in Fowler et al.| ( |2010| . Here, for every single time- 
stream sample cut, we fit for its (unknown) amplit ude as part 
of the map solution, as suggested by, for example, |Patanchon] 
et al. ( 2008| l. In addition to cuts from cosmic rays and the 
like, we explicitly cut the first and last second of each TOD 
to give the mapping algorithm the freedom to match the TOD 
beginning and end as smoothly as possible. Since these data 
are already highly downweighted by the time-domain taper, 
the additional amount of data loss is negligible. 

11.4. Map Solution 

With these pre-processing steps, form of the noise, and 
form of the solution specified, we must then actually solve 
for the map. We use the classic Preconditioned Conjugate 
Gradient algorithm, using the hit-count map as a Jacobi (di- 
agonal) preconditioner for the sky map part of the solution. 
To recover the scales around I ~ 200 and below, we need a 
few hundred PCG iterations. The map making is quite CPU- 
intensive, with each iteration taking about a minute on 1600 
2.53GHz Nehalem cores on the Scinet GPC, or a full run to 
1000 iterations taking about a day of wall clock time, or 5 
years of CPU time. 

We note that care must be taken when estimating the noise 
to make sure it does not bias the map. Consider the follow- 
ing thought experiment: Two detectors are scanned across 
a source, and their weights are estimated from the internal 
scatters of their time streams. In general, noise in the de- 
tectors will lead to one of them observing a lower flux from 
the source. If the detector weights are then measured and 
the source is not removed from the time streams, the detec- 
tor which happened to measure the lower flux will usually 
have a lower variance in its data, and will therefore on aver- 
age receive more weight than the other detector. When the 
detectors are combined, the more heavily weighted low-flux 
measurement will lead to the source flux being systematically 
biased low. The magnitude of the bias is set by the SNR 
in the individual time streams, not the final SNR. The bias 
can be mitigated if an estimate of the signal is removed from 
the time streams before noise estimation. Since the CMB is 
highly subdominant to the noise in ACT time streams, this 
effect is small, but simulations show that it is not negligible. 
So, first we make an initial map using the full dataset and the 
noise measured directly from the time streams, then we sub- 
tract that map from the time streams to reduce the sky signal 
in them before estimating the noise for a second time, and fi- 
nally we use this improved estimate of the noise to make the 
final maps. 

ACT must take a bit more care than all-sky experiments be- 
cause of the large change in sensitivity across the map. In 
particular, near the map edges, only a small amount of data 
contributes to the map, and so removing the raw map would 
lead to an artificial reduction in the noise of the data near 
the edges, which would lead to an artificial increase in their 
weight. To prevent this from happening, the starter maps used 
are first filtered and apodized. The apodization is generated 
by first rescaling the hit count map to the 0-1.0 range, then 
setting all values higher than a threshold to unity, and finally 
smoothing the resulting weight map with a Gaussian window. 
Furthermore, from the starter data map we filter out scales 
larger than I = 300 where atmospheric noise is very large, and 
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scales smaller than I = 3000 where the map is highly noise- 
dominated and simulations show that the bias effect is truly 
negligible. Finally, we generate the starter map by multiply- 
ing the the filtered data map with the apodization window. To 
test this procedure, we carry out the two-step procedure on 
both the real data, and the real data with a simulated map in- 
jected into it. We compare the difference of these two sets of 
maps to the original simulation, and show the resulting trans- 
fer function in Figure 10 At the low-pass scale of I = 300, the 
maps are unbiased to 03%, and rapidly become unbiased to 
better than a part in 10 3 on smaller scales. 

1 1 .5. Map Analysis and Results 
The resulting maps cover and area of 845.6 deg 2 on the sky, 
ranging between 20 h 43 m and 7 h 53 m in right ascension and 
-48? 1 and -57?2 in declination. A total of five maps were 
made: one using the full dataset, plus four partial versions us- 
ing independent subsets of the data for noise estimation pur- 
poses. Each map is accompanied with a hit map with the 
number of data samples that fell on every pixel of the cor- 
responding map. 

Thes e maps are a n overall im proveme nt ov e r those used in 
|Fowler et a"L](|2010|i;[M arriage et al.|(|201 la|b|l;|Dunkley et al. 
(1 201 1 J ) ; |Da s e t al . | ( f 20 1 1 b[ l ; |Hajian et al.| ( |201 1| ; |Sehgal et al. 
(201 1 ), marked particularly by lower noise on angular scales 
between 500 < I < 2000, reducing the variance by up to a fac- 
tor of 2 or 3. However, the original maps were already nearly 
sample-limited on these scales, so errors on the power spectra 
are only modestly reduced. These changes are due to slightly 
different data selection cuts, improvement in the noise treat- 
ment during the mapping process and improved per-detector 
calibration. 



These sensitivities were computed using the number of sam- 
ples per unit area and the noise equivalent temperature given 
in Table [3] We run for 2000 PCG iterations on this map to 
ensure that it is truly converged, though as mentioned we find 
that all science scales are converged in far fewer PCG itera- 



tions. 



Figure 12 shows the power spectrum computed for this 



map, c ompared to previ ous versions a lready published in 
Fowler et al. ( 2010| > and Das et al. ( 201 lb| l. The uncertainty 
in the spectrum reveals an improvement in the quality of the 
map. 

The parameters obtained using a fit to the spectrum released 
with this pa per are slightly s hifted with respect to the values 
presented in Dunkley et al. (2011 1. The largest shifts in the 



primary parameters were for n s and VLi,h 2 at Act = 0.3 and 
with Act = 0.7, although the change in the inferred dark en- 
ergy density changed by a negligible 0.03 a. The other param- 
eters changed by less than 0.07 a. The secondary parameters 
were affected by the new mapping procedure, decreasing the 
Poisson point source level from 13.7/iK 2 to 12.5 fiK 2 . The 
upper 95% confidence limits of the correlated point source 
amplitude and SZ amplitude are instead unchanged. 

Also, a cross-correlation analysis to the BLAST maps 
( |http : //blastexperiment . inf o/) show correla- 
tions detected at a 25 a level, implying a detection of emission 
by radio and du sty star-form ing galaxies (DSFGs) at high I, 
as described in Hajian et al. ( 2012) . 

Finally, the flux of sources as estimated in the data release 
is approxim ately 5% higher than those reported in Marriage 
|et al.| ( |201 la| l, due to the new source treatment in the map 
making and a change in the beam solid angle. The best repre- 
sentation of the beam shape, as well as the window function 
for power spectrum analysis, are released together with the 
maps. 
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FIG. 10. — Transfer function of the map making process. A simulated map 
(in) was added to the time streams and mapped together with the real data 
(sim-inject map). Then the real data map (previously computed) was sub- 
tracted from the sim-inject map to produce the output map (out). The transfer 
function was computed as the power spectrum of the output map divided by 
the cross-spectrum of the output map and the simulated map. The result is es- 
sentially unity for all scales above our noise-estimation input map filter scale 
of 300 (see text), with a slight boosting of large scales that reaches 2% at 
multipoles of 200. 



Figure 1 1 shows the final 2008 southern map, overlaid by 
contours of iso-sensitivity (35, 50 and 65 /iK-arcminute), as 
different regions of the map have different integration times. 



12. CONCLUSIONS 

We have fully characterized the data from the 2008 season 
of observations from ACT, including data selection, calibra- 
tion, pointing, random and systematic noise and atmosphere 
contamination. 

Observations done in 2008 yielded a total of nearly 
1260 h of CMB survey data at 148 GHz, distributed between 
two observation stripes centered at declinations near 0° and 
-53°, with approximately 850 and 280 square degrees respec- 
tively. After data selection, the observing time was reduced to 
nearly 816 h. Out of a potential 1024 live detectors, the num- 
ber of effective detectors was on average 593. Combining 
observing time and detector performance, it yields a system 
efficiency of 38%. Including the ratio of observing time over 
calendar time, the overall system efficiency was 15%. 

The uncorrelated noise in the data, which excludes the at- 
mospheric fluctuations, is dominated by detector noise. Be- 
tween 5 Hz and 20 Hz, we found a total NET of 32 /iK-^/s at 
148 GHz. In other words, given 1 h of observation per square 
degree, or 1 s of integration per square arcminute, the noise 
should be about 32 /iK arcminutes. 

The correlated noise is dominated at low frequencies by 
the atmosphere and thermal drift, while at higher frequen- 
cies electromagnetic pickup and mechanical vibrations are the 
largest sources of correlated noise. The former are character- 
ized in frequency by negative power laws, equaling the detec- 
tor white noise spectra at approximately 2 Hz. The latter show 
up mostly as narrow-band signals emerging from the detector 
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FIG. 11. — Final map of the southern region observed at 148 GHz in 2008 obtained after 2000 PCG iterations. The contours in white join regions with the 
same sensitivity, namely 35, 50 and 65 /^K-arcminute, with increasing sensitivity towards the center of the maps. The total area enclosed by the contours is 
220 deg 2 , 420 deg 2 and 530 deg 2 respectively. The region between 18 A 54 m and 2\ h l5 m in right ascension is not shown because it was sparsely observed and has 
low sensitivity. Most of the analysis has been done in the deeper region between roughly 0* and 7* in R.A.. The map has been high-pass filtered for clarity, 
depressing modes larger than i = 300. The CMB anisotropics can be seen by eye. There is some evidence of large scale systematic noise, especially near the 
edges, observed as long horizontal features, which are mostly related to scan-synchronous systematics (such as ground spillover). Despite the size of the plot, it 
is still possible to see some bright point sources and SZ clusters. 
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FIG. 12. — Power spectrum of the 148 GHz map obtained after 2000 PCG iterations compared to the power spectra in [Fowler et al.|(2010) andJDas et al.| 
(201 lb) . The thick orange curve shows the best-fit model including the CMB secondaries and point-source contribution taken from Uunkl e~et al.|pOl 1^ . As 
revealed by the error bars, the power spectrum derived from this map represents an improvement from previous maps. 




FIG. 13. — Side by side comparison between the ACT map (season 2008) and the SPT map for the same region of the sky. The left panel shows the ACT map 
high-pass filtered with a cos 2 £-like filter that goes fro m to 1 for 100 <£ < 300, and the center and right panels show the ACT and SPT maps respectively under 
the same high-pass filter used in the SPT data release I Schaffer et al. 201 1 1. Agreement between the CMB features in the two maps is clear by eye. 
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noise. They can all be modeled as correlated modes added to 
each detector TOD. Once identified, they can be de- weighted 
in the map solution, minimizing the errors in the map, while 
keeping the solution unbiased. 

The data processing steps for map making can be summa- 
rized as follows: The median signal level of each detector is 
removed together with a single array-wide slope across the 
15-minute file, and the inverse anti-alias filter and time con- 
stant deconvolution are applied (see ^5] and 5|9ji. Then a pre- 
calculated data selection, calibration and pointing solution are 
applied (see |T0] ^7]and ^8j. After this, 7 pre-computed dark 
modes are de-projected from each detector TOD (see 5 1 1 and 
Appendix [A|. Before map making, the expected signaTfrom 
the bright point sources is subtracted from the time streams. 
Then the maps are made by minimizing \ 2 over a noise model 
which includes correlated modes, cuts (missing time streams) 
and frequency dependence, as described in JIT] This is done 
twice, with an estimate of the map removed from the data be- 
fore the noise estimations are repeated in the second mapping 
pass. Finally, the flux from point sources previously removed 
is re-added to the maps. The resulting maps cover 845.6 deg 2 
on the sky and are consistent with WMAP at angular scales 
measured in common. 

The same area of the sky covered by these maps was ob- 
served by the South Pole Telesco pe (SPT) team, wh o have 
recently made their data public (Schaf fer et al.| |201 1 ). Fig- 
ure [13] provides a side by side comparison between the ACT 
map and the SPT map for the same region of the sky and us- 
ing the same filtering used in the SPT release. The correla- 
tion between CMB features in the two maps is clear to the 
eye, supporting the excellent quality of both measurements. 
It is also clear that the ACT map is noisier than the SPT map, 
which will improve in future map releases when data from the 
following observing seasons (from years 2009 and 2010) are 
included. 

Sky maps for the 148 GHz ACT southern observations from 
2008, described in this paper, are available through NASA's 
Legacy Archive for Microwave Background Data Analysis 



(LAMBDA), where a variety of ACT analysis software, data 
products, and model templates are also available. Future data 
releases will include ACT observations at higher frequencies 
and subsequent observing seasons, as well as sky coverage in 
ACT's equatorial stripe which overlaps numerous other ob- 
serving programs. 
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APPENDIX 

MODE SELECTION AND REMOVAL 

Contaminant signals, like the atmosphere and systematic signals, produce correlations between detector TODs. Removing this 
"correlated" noise is important for both map making and data characterization, so we devote this appendix to explain how this is 
done in more detail. 

A contaminant signal can be modeled as a time stream superimposed on the TODs from a set of individual detectors (common 
mode). We call this time stream a "correlated mode." One way to estimate it is using an appropriate linear combination of 
detector TODs. For instance, the common mode, defined as a linear combination with equal weight of all detector TODs, is a 
good estimator of the atmosphere signal. Organizing all the detector TODs into columns of an n samp \ e by ridetector matrix A, the 
correlated mode m can be expressed as 



m 



(Al) 



where v is a vector of unit magnitude representing the nd elec tor coefficients of the linear combination, and s normalizes the mode 
so it also has unit magnitude. We henceforth define the correlated modes such that they are always normalized. 
Once a correlated mode has been identified, it can be fitted and subtracted from the data: 



A': 



I-mm' A 



:A(T-vv r A r As- 2 ) 



(A2) 



We can also construct correlated modes m in other ways besi des t he linear combinations in Eq. Al (for example, from a 
thermometry time stream), and in this case the first equality in Eq. A2 must be used directly, rather than the second equality. 
One can identify important correlated modes using the singular value decomposition (SVD) of the matrix A: 



A = USV r , 



(A3) 
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where the columns of U can be identified as normalized corr elate d modes m, S is a diagonal matrix containing the singular values 
of A (identified as the normalization factor, s, in E quation \Al\ , and the columns of V are the eigenvectors of the covariance 
matrix A r A (identified as vector \ in Equation [aT]). The SVD mode selection can be tailored to specific signals by first finding 
the common mode from subsets of detector TODs, and then applying the previous method to find the strongest correlated modes 
out of this reduced set of modes. For example, to find modes that are more likely to correlate in rows, we first find the common 
mode from all detector TODs in each row (32 row common modes in total), and stack them together as columns of a matrix A rmv . 
The strongest modes are then readily found using SVD. This method is useful to identify row- and column-correlated modes, and 
sub-array-scale modes from the atmosphere. For the latter, we divide the array into 16 square blocks, finding the common mode 
in each one before applying the SVD. We call the set of modes found this way a "multi-common mode." 

An important consideration when using modes obtained from linear combinations of live detectors is that they are also corre- 
lated to the sky signal. For example, naively removing the common mode effectively filters the sky, depressing scales larger than 
the array size (24' or £ < 450). Naturally, the effect is stronger for the multi-common mode, which filters scales larger than a 
fourth of the array (6' or I < 1800). For this reason, the multi-common mode is used to calculate the drift error, but it cannot be 
naively removed for making maps, nor can any other mode obtained from live detectors. 



In Section 1 1 the appropriate way of projecting modes out of the data is discussed in the context of map making. 



Dark Mode Removal 

The dark detectors share the same readout as the live detectors, so they can capture systematics like thermal drift, electromag- 
netic pickup and magnetic pickup. Moreover, as they are not coupled to the optical signal, correlated modes obtained from dark 
detectors (dark modes) can be safely removed as a preprocessing step before making maps. For example, the thermal drift is well 
represented by the low-pass-filtered common mode of the dark detectors. 

For row- and column-correlated electromagnetic signals, the ability to identify correlated modes depends on whether the 
desired row or column of the array contains a working dark detector. The array design uses one dark detector per column, 
always in the same row. This is ideal for identifying column-correlated modes, but not for row-correlated modes. For the latter 
we must rely on the broken live detectors with a properly working readout circuit, as described in ^ flO] Nevertheless, the row- 
and column-correlated modes can still be identified by using SVD analysis over all dark detectors. Note that it is useful to first 
remove the slow thermal drift from the dark detectors before trying to find these other higher frequency signals. Figure [7] shows 
an example of correlation matrices before and after having removed the dark modes. Note that both column and row correlations 
are significantly suppressed after the removal. 
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